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Generation of Bell, W and GHZ states via exceptional points in non-Hermitian 

qnantum spin systems 
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We study quantum phase transitions in non-Hermitian XY and transverse-field Ising spin chains, 
in which the non-Hermiticity arises from the imaginary magnetic field. Analytical and numerical 
results show that at exceptional points, coalescing eigenstates in these models close to W, distant 
Bell and GHZ states, which can be steady states in dynamical preparation scheme proposed by T. 
D. Lee et. al. (Phys. Rev. Lett. 113, 250401 (2014)). Selecting proper initial states, numerical 
simulations demonstrate the time evolution process to the target states with high fidelity. 
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I. INTRODUCTION 

Quantum phase transition can occur in a finite non- 
Hermitian system, associating parity-time {'P'T) reversal 
or other type of symmetry breaking. At the transition 
point as referred to exceptional point (EP), a pair of 
eigen states coalesces into a single state. Many finite¬ 
sized discrete systems have been investigated, including 
tight-binding models, quantum spin chains, and complex 
crystal. 

These features are different from that of quantum 
phase transition in infinite Hermitian system. Recently, 
critical behavior of non-Hermitian system has been em¬ 
ployed to generate entangled states in a dynamical pro¬ 
cess and the corresponding experimental protocol is also 
proposed [l|, [3 . According to the non-Hermitian quan¬ 
tum theory [SHlOj, a pseudo-Hermitian system has real 
eigenvalues or conjugate pair complex eigenvalues. Con¬ 
sidering the simplest case, there is only a single pair of 
eigenstates breaking the symmetry of the Hamiltonian, 
with conjugate complex eigenvalues. A seed state is an 
initial state consisting of various eigenstates with eigen¬ 
values with zero, positive and negative imaginary parts, 
respectively. As time evolution, the amplitude of the 
state with positive imaginary part in its eigenvalues will 
increase exponentially and suppress that of other compo¬ 
nents. The target is the final steady state and expected to 
have peculiar features for quantum computation process¬ 
ing and other applications. It is important to construct 
a simple Hamiltonian which is suitable for experimen¬ 
tal implementation: to prepare desirable quantum states 
with high fidelity. 

In quantum information science, it is a crucial problem 
to develop techniques for generating entanglement among 
stationary qubits, which plays a central role in applica¬ 
tions |lll - ll3| . Bell states are specific maximally entan¬ 
gled quantum states of two qubits. For many-qubit sys¬ 
tem, there are two typical multipartite entangled states, 
Greenberger-Horne-Zeilinger (GHZ) and W states, which 
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are usually referred to as maximal entanglement. Multi¬ 
partite entanglement has been recognized as a powerful 
resource in quantum information processing and commu¬ 
nication. Numerous protocols for the preparation of such 
states have been proposed [3®. 

In this paper, we consider whether it is possible to use 
non-Hermitian systems to generate a W, distant Bell and 
GHZ states via the dynamical process near EPs. We in¬ 
troduce a non-Hermitian XY and a transverse-field Ising 
spin chains to demonstrate the schemes. Numerical simu¬ 
lations show that the target states can be obtained with 
high fidelity by the time evolutions of selecting proper 
initial states. 

The remainder of this paper is organized as follows. 
In Sec. im we present a non-Hermitian XY spin model 
and solutions. Secs. nil 113 and El are devoted to the 
schemes of preparing W, Bell and GHZ states, respec¬ 
tively. Finally, we present a summary and discussion in 
Sec. ED 

II. XY SPIN CHAIN 

We consider a non-Hermitian XY spin model 

N-l 

i?chain = 2 X + H-C. 

{V+ i-i)al + {V ( 1 ) 

on an A-site chain, where af {a = x, y, z) is Pauli matrix. 
In the case of 7 = 0, it is reduced to a Hermitian model 
with V symmetry. Here the parity operator V is given by 
= (jf with T = {N+l — l). In the case of nonzero 
7 , the V symmetry is broken, but VY is still symmetric, 
where T is a time reversal operator 'TiT~^ = —i. 

We note that 

[•^z,chain] = 0, (2) 

where Ja = ^ total spin operator. This means 

that iLchain Can be diagonalized in each invariant sub¬ 
space.In this paper, we only concern the issue in the sub¬ 
space with Jz = N — 1 and N =even. In this invariant 
subspace, the wave function has the form 
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FIG. 1. (Color online) Plots of the fidelity / (t) for preparing W state, as a function of time for the systems with = 6 (a) 
and 8 (b). The times are dimensionless and in units of 1. We see that the fidelities converge to constants in exponential manner 
and the converging hdelities become higher as 7 closes to 1, while the converging times get longer. The obtained results are 
not sensitive to the size N, which is quite different from the situations for the productions of Bell and GHZ states. 


N 

= (3) 

where |JJ-) is a saturated ferromagnetic state |J|) = 
Uhli)- Then we get an equivalent Hamiltonian 


N-l 

Heq = 'y ] \l) (Z + 1| + H.C. 

1=1 

+ {V + i^)\l){l\ + {V-i^)\N){N\, (4) 

where the position state at Zth site is |Z) = |JJ.). The 

eigen problem of the equivalent Hamiltonian is given in 
Appendix. In the following, we will discuss the schemes 
for the preparation of W and Bell states based on the 
Hamiltonian Heq. 


III. W STATE 

In the situation V = 0, the Hamiltonian Heq is reduced 
to 


N-l 

Hw=Y1 IO(^ + l|+H.c.+l 7 |l)(l|-l 7 |iV)(iV|. (5) 
1=1 

The exact solution in the Appendix suggests us to con¬ 
sider the state 

N N 

= ( 6 ) 

which represents a single-magnon spin wave with wave 
vector tt/2. It is a W state under a local transforma¬ 
tion {—iy cr/" —>■ a-y, which does not reduce its properties 


in quantum information processing. A straightforward 
derivation shows that the state |W) is an eigenstate of 
the Hamiltonian i?w at 7 = 7 c = 1, i.e., 

Hw ( 7 c) |W) = 0 |W). (7) 

Then, we will show that |W) is a special eigenstate of 
Hwilc)- For the corresponding conjugate Hamiltonian 
( 7 c), we have 

K (7c) |W) = 0 IW), ( 8 ) 

where 

1 ^ 

= ( 9 ) 

It is easy to find that 

(w |W) = 0, (10) 


which indicates that Hw Fas an EP at 7 c and the W 
state is the coalescent state at the transition point. In 
the Appendix, this result is confirmed by an exact Bethe 
Ansatz analysis. 

We investigate the scheme of selecting the |W) state by 
a dynamic process. From the Appendix or the previous 
work [m, we find that the complex conjugate pair of 
energies is ±f 2 JsinhK for small k, where real number k 
obeys the equation 

7 ^ sinh[(A — 1 ) k] = sinh[(A -|- 1 ) k]. (11 ) 


We note that the value of 7 determines the gap between 
the complex conjugate pair of energies, or the converg¬ 
ing time. The initial state is taken as |'0 (0)) = |1), the 
evolved state lip {t)) is expected to close the target state 
for sufficient long time. We employ the fidelity 


fit) 


(W|^(t)) , 


( 12 ) 
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to characterize the efficiency of the scheme. Here 




is the Dirac normalized state of \iIj (<)) to reduce the in¬ 
creasing norm of \^p (t)). In the limit case of 7 —>■ 1, we 
will have / (t) —>■ 1 as t —>■ 00 . For finite 7 , the time 
evolution of the state is computed by numerical diago- 
nalization in the broken symmetric region. In order to 
quantitatively evaluate the fidelity and demonstrate the 
proposed scheme, we simulate the dynamic processes of 
the W state preparation. To illustrate the process, we 
plot the fidelities as functions of time for systems with 
N = 6 and 8 in Fig. [TJ It shows that the fidelities con¬ 
verges to a steady value exponentially fast. Smaller 7 
(approaches to 1 ) can enhance the fidelity, while the con¬ 
verging time becomes longer. Moreover, we find that the 
converging times for two cases are not so sensitive to the 
size iV, which is quite different from that in following two 
schemes for preparing distant Bell and GHZ states. This 
is because of the fact that the phase boundary is always 
at 7 = 1 for any even N. Then such a scheme is more 
efficient for a W-state production. 


IV. BELL STATE 

In the situation \V\ > 2, the exact solution in the 
Appendix shows that two bound states are formed, in 
which the probability mainly distributes around two end¬ 
ing sites. The phase diagram has been obtained as Eq. 
(IA.I 6 I) in the Appendix, which is the base of the scheme 
for preparing Bell state. According to the Bethe Ansatz 
result, there is a conjugate complex pair of energy levels 
in the broken VT symmetric region. The magnitude of 
the imaginary part of the eigen energy |Ime| is also an 
indicator of the phase boundary and determines the con¬ 
verging speed of the scheme. For illustrating this point, 
we plot |Ime| as a function of V and 7 for the systems 
with iV = 6 , 8 and 10 in Fig. [21 The corresponding 
exact boundary from Eq. (IA.I 6 I) is plotted as well. We 
find that they accord with each other and the boundary 
appears as a linear line with V-dependent slope in the 
logarithm scales. We will see that the profile of the phase 
diagram directly determines the efficiency of the scheme 
in the following investigation. 

In order to understand a clear physical picture of the 
exact solution, we use the perturbation method to sim¬ 
plify the Hamiltonian Hgq in large V limit. Although 
the perturbation theory for non-Hermitian Hamiltonian 
has not been well established, the following result will 
show that the corresponding approximation is technically 
sound by the comparison with the exact solution. We 
rewrite the Hamiltonian Hgq in the form 

Heq = Ho + H', (13) 

N-2 

Ho=Y^ \l) (;-bI|+H.c. 

1^2 

+ iV + zj)\l){l\ + iV-ij)\N){N\, (14) 

iI' = |l)(2| + |V-I)(iV|+H.c., (15) 


where the eigen states of Hq can be easily obtained as 

{|1), \N), Ef=-/ ^ e [2,iV- 1]} 

with corresponding energy {V + i"/, V — i'f, 2cos[^^^(^]; 
n G [2,iV — 1]}. This set of eigen states has a special 
feature that they can construct a complete set under the 
Dirac inner product, even Hq is a non-Hermitian Hamil¬ 
tonian. Then the effective Hamiltonian for two bound 
states can be obtained as 


Heff = Aeff |1) {N\ + H.C. + {V + Eeff + ^7) |1) (1| 

+ iV + V,s-ij)\N)(N\, (16) 

in the case of \V\ :$> 1, the model above is a simple two- 
site model and easily solvable. Here the effective poten¬ 
tial is 


- AT _ 1 E 


o N-1 . 2 

2 sin < 


n—2 


1 


V — 2 cos (/)„ V 


(17) 


and the effective coupling is 

2 sin 4>n sin[( A — 2 ) 


N — 1 ^ V — 2 cos ( 

n—2 

n 

V2’ 


(18) 


where parameters O, (jzn and 9 are N dependent functions 


^ cos [{N - 4) 7 r/ 2 ] sin[(A - 4) (A - 2) 9] 
(A - 1) sin[(A - 4) 9] 

(-l)^/%in[(A- 2) A 6 »] 

(A - 1) sin(A 6 l) ’ 

(j)n =2{n- 1)9 

9 = ^ 

2(A- !)■ 


(19) 

( 20 ) 
( 21 ) 


The eigen states of Aes are 


(*7±7A2g-72)|l) + Aeff|A), (22) 


with eigenvalues: ± A^fj — 7 ^ + V + VeS- At the EP, 
Agfj = 7 ^, the coalescent state is 


z7e|l) + Aeff|A), (23) 

with energy 

Sc = V + Ves ~ ^ A y j (24) 

which is in agreement with the approximate expression 
Eq. (IA.19I) in the Appendix. Then the boundary has the 
form 


In I 7 I -I- 2 In I V| = In |H| 


(25) 
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FIG. 2. (Color online) Phase diagram of non-Hermitian Hamiltonian in Eq. (jd)). The color contour map represents the 
magnitude of the imaginary part of conjugate-pair energy levels for = 6 , 8 and 10, obtained by exact diagonalization. The 
white area indicates the region where the spectrum is entirely real. The dot line is the plot of function in Eq. (IA.16I) . indicating 
the exact phase boundary. We see that all three boundaries are in line shape with different slopes in the logarithm scales. This 
property can be explained by perturbation approximation. 



FIG. 3. (Color online) Plots of the fidelity / (t) for preparing Bell state, as a function of time for the systems with = 6 (a) 
and 8 (b). The times are dimensionless and in units of 10® and 10"^, respectively. We see that the converging fidelities approach 
to 1, getting higher as increasing of V and longer time. The time scale of (b) is over ten times longer than that of (a), which 
indicates the difficulty of preparing long-distance Bell state. 


in the logarithm scales, indicating a linear phase bound¬ 
ary with a fixed slope. This is qualitatively in agreement 
with the numerical results in Fig. [^obtained by the ex¬ 
act solution, where the slopes of the boundary are N 
dependent. 

Based on the phase boundary, one can prepare the tar¬ 
get state in the vicinity of the EPs via dynamic process. 
The target state is a Bell state, expressed as 

|Bell) = ^(|l)-*|iV)). (26) 

The initial state is taken as |'0 (0)) = U); the evolved 
state \ip {t)) is expected to close the target state for suf¬ 
ficient long time. We employ the fidelity 


fit) 


(Belli ^(t)) , 


(27) 


to characterize the efficiency of the scheme. Here 


f’it) 


is the Dirac normalized state of \'tp (t)) to reduce the in¬ 


creasing norm of {ip {t)). The time evolution of the state 
is computed by numerical diagonalization. For given N 
and V, we numerically search an optimal 7 to obtain 
higher fidelity in the broken symmetric region. In order 
to quantitatively evaluate the fidelity and demonstrate 
the proposed scheme, we simulate the dynamic processes 
of the quantum state preparation. To illustrate the pro¬ 
cess, we plot the fidelities as functions of time for systems 
with N = 6 and 8 in Fig. [31 It shows that the fidelities 
converges to a steady value exponentially fast. Larger 
\V\ corresponds to smaller optimal 7 , leading to higher 
fidelity, but longer converging time. We also find that the 
converging times for two cases are sensitive to the size N. 
These accord with the phase diagrams in Fig. linear 
boundary indicates that larger ln|F| matches smaller Iny 
and slight change of slopes between ln|F| and Iny results 
in drastic change of the converging times. 
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FIG. 4. (Color online) Phase diagram of non-Hermitian Hamiltonian in Eq. (I29II . The color contour map represents the 
magnitude of the imaginary part of conjugate-pair energy levels for N = 6 and 8, obtained by exact diagonalization. The white 
area indicates the region where the spectrum is entirely real. We see that two boundaries have similar shapes but with a shift. 



FIG. 5. (Color online) Plots of the fidelity / (t) for preparing GHZ state, as a function of time for the systems with = 6 (a) 
and 8 (b). The times are dimensionless and in units of 10® and 10“*, respectively. We see that the converging fidelities approach 
to 1, getting higher as decreasing of A and longer time. The time scale of (b) is over ten times longer than that of (a), which 
indicates the difficulty of preparing long-distance GHZ state. 


V. GHZ STATE 

The above conclusion provides a way to prepare a su¬ 
perposition of two distant position states. Such a scheme 
can be extended to prepare the GHZ state which has the 
form 

N 

|GHZ) = |^)+n'^,+ l^)- ( 28 ) 

States in) and Ilzli 1^) can be regarded as two end 
position states, which are connected by A^-step opera¬ 
tions of operator ^i- This opens a probability to 

select the GHZ state as a steady state near the EP. We 
consider a simple and practical model, which is a non- 
Hermitian Ising model, described by the Hamiltonian 


N N N 

Hghz = -JY1 + *7 E ^ E ■ (29) 

1=1 1=1 1=1 

It is a standard transverse-field Ising model at 7 = 0, 
which can be exactly solved and has been extensively 
studied in a variety of areas. Recently, theoretical stud¬ 
ies of several types of quantum Ising models were ex¬ 
tended to the non-Hermitian reg ime and some peculiar 
properties were observed [32 - 1^ . In the case of J = 0, 
this model is reduced to non-interacting spin-I /2 par¬ 
ticles with complex magnetic field, which has full real 
spectrum when ^ 7 ^ [ 3 ^. We assume that the phase 
transition can occur in the case of nonzero 7 and J. Since 
this model is not solvable, we perform numerical simula¬ 
tion by exact diagonalization. 
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Similarly as last section, we still employ the magnitude 
of the imaginary part of the eigen energy |Im£| as an 
indicator to characterize the phase boundary. Taking 
J = 1, we plot |Ime| as a function of V and 7 for the 
systems with TV = 6 and 8 in Fig. |4l We find that the 
phase boundaries of the two cases have the similar profile 
but with a shift. This will be reflected on the speed of 
the fidelity convergence. 

For a GFIZ state preparation, the initial state is taken 
as IV'(0)) = | 1 )) the evolved state |V'(t)) is expected to 
close the target state for sufficient long time. We employ 
the fidelity 


f(t)= (GHZ|^(t)) 


(30) 


to characterize the efficiency of the scheme. Here 


V’(t) 


is the Dirac normalized state of |V’ (i)} to reduce the in¬ 
creasing norm of {ip {t)). The time evolution of the state 
is computed by numerical diagonalization. For given N 
and V, we numerically search an optimal 7 to obtain 
higher fidelity in the broken symmetric region. In order 
to quantitatively evaluate the fidelity and demonstrate 
the proposed scheme, we simulate the dynamic processes 
of the quantum state preparation. To illustrate the pro¬ 
cess, we plot the fidelities as functions of time for systems 
with N = 6 and 8 in Fig. [SJ The obtained results are 
similar to the case of Bell-state production at last. 


VI. SUMMARY 


and the EPs in the cases of U = 0 and \V\ > 2. 


1. V = 0 case 


The Bethe Ansatz wave function is in the form 

N 

\k) = Y, |j), (A.2) 

f=i 

where k is real number, indicating a scattering state. The 
Schrodinger equation H \ k) = Sk\k) can be written as 


M 


Ak 

Bk 


where the matrix 
M = 


v+e^^ + 


= 0 , 


v+e -I- e 


(A.3) 


^^^ikN _j_ gife(iV-l) ^^^-ikN _j_ ^-ik(N-l) 

v± = ±ij-ek, (A.4) 

and the real spectrum 

efc= 2 cosfc. (A.5) 

The existence of solution requires 

det |M| = 0, (A. 6 ) 

which leads to the equation 


In summary, we presented schemes to generate W, dis¬ 
tant Bell and GHZ states by exploiting the quantum 
phase transitions in non-Hermitian XY and transverse- 
field Ising spin chains. The phase diagrams for such two 
models are obtained analytically and numerically, which 
is crucial for the practical realization of the scheme. Nu¬ 
merical simulations on the dynamics process for state 
preparation show that the evolved states close to tar¬ 
get states in an exponential manner over time. Gom- 
paring the dynamical preparation of quantum state via 
Hermitian system, where the acquired state only emerges 
within a short time window, this scheme can provide the 
steady final state. A shortcoming of the scheme is that 
the production period for Bell and GHZ states increases 
rapidly as cluster size grows. However, this scheme is 
more efficient for a W-state production. 

Appendix: Exact solution of the Heq 

In this appendix, we present the exact results for the 
solutions of following model 

AT-l 

iVeq = 10 + 1| + H.C. 

1=1 

+ (U + *7)|l)(l| + (H-*7)|TV)(iV|, (A.l) 


F (fc) = sin [k {N + 1)] -|- 7 ^ sin [k {N — 1)] = 0. (A.7) 

The EP kc can be determined by equation 

F{k,) = ^F{k,) = 0. (A. 8 ) 

We obtain fee = 7 r /2 at 7 = 1 (see ref. [3l|). 

2. |V| > 2 case 

In this situation, we are interested in bound states. 
The corresponding Bethe Ansatz wave function is in the 
form 

N 

|«) = ^(a„e«^+/3,e-«^)|j), (A.9) 

where k is a real number. By the similar procedure, we 
reach the equation 

F (kc) =-^F (k,) = 0, (A.IO) 

which determines the location of EPs at energy 
= 2cosh«;c, 


(A.ll) 
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where function 

F (k) = sinh [(iV + 1) k] — 2V sinh (Nk) 

+ + 72 ) sinh [(TV - 1) k] . (A.12) 

From Eq. (lA.lOl) . we have 

[Nt]+ + ? 7 _] coshKc — 2NV 
{Nri_ + 77 +) sinh Kc 
rj- sinh 


C IS 


= — tanh (A^Kc) , (A.13) 


77 + coshKc — 2V 

where 

? 7 ± = 1±E^±7^. (A.14) 

The bound state EPs require 

|e«ol>2|1^|. (A.15) 


c = coshK = F [1 


4A^ (?7+ - 1) + 1+1- + 4A^F^) 


F = 


y2 {2Nr)++r]_Y 


V{2Nr]+ + rj_) 


4iV(r7+-l) 

In the case of \V\ 1, we have 

V 1 


(A. 18 ) 


(A.19) 


which gives the approximate energy expression = 
2 coshK Ri E + A. 


Such solutions exist when parameters V and 7 satisfy 
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